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I. INTRODUCTION 



The Schwinger model is QED in 1+1 space-time dimensions 0]. Although it is a toy model, it shares many 
interesting physical properties with more realistic theories such as QCD or the electroweak theory. It is perhaps 
the simplest example in which gauge invariance does not necessarily imply a massless gauge boson, analogously to 
the Higgs phenomenon. Other interesting properties of the model are dynamical mass generation, chiral symmetry 
breaking and confinement. The model with massless fermions was shown to be exactly solvable in vacuum (that is, 
without thermal effects) long time ago [0j2|. It is equivalent to a theory describing a free boson with mass ejypn 
(bosonization), which is physically a fermion-antifermion bound state (confinement). The finite mass implies a finite 
correlation length, which physically corresponds to charge screening, long range forces being absent. On the other 
hand, chiral symmetry is broken through the chiral anomaly H rather than spontaneously, since Coleman's theorem 
Q prevents any continuous symmetry to be spontaneously broken in two dimensions. When a mass parameter for 
the fermions is included, the model is no longer solvable but it is still possible to analyse exactly some of the above 
properties, like fermion confinement || . 

In the last ten years, there has been a renewed interest in the study of the Schwinger model including statistical 
mechanics features. The model was also solved and, in particular, the two point correlation functions and the partition 
function were obtained at finite temperature T and zero chemical potential fx in || when the Dirac operator has no 
zero modes. Its thermodynamics can be expressed in terms of those of a free boson of mass e/^/ir and free fermions, 
i.e., bosonization also takes place at finite temperature. A more complete study of the Schwinger model on a torus, 
which naturally incorporates temperature effects, still at fx = 0, has been performed in @,^]. More specificallyin Q] 
the model was treated with an arbitrary number of zero modes and the two-point fermion correlation function was 
calculated, whereas in || higher correlation functions were obtained. In the correlation functions have also been 
studied for nonzero fx. In a recent work |l0| , the problem of charge screening at finite temperature (with fi = 0) in 
the Schwinger model is analysed, in connection with the spontaneous breaking of the discrete Z symmetry, which 
corresponds to the freedom of choosing gauge fields in the Euclidean time direction with any winding number around 
S . These nontrivial gauge transformations will play an essential role in our analysis with a nonzero chemical potential. 

On the other hand, we remind that the Thirring model, which describes massless fermions in 1+1 dimensions with a 
quartic self-interaction, can also be explicitly solved in vacuum (T = fx = 0) []Tl 12 1. We also recall that the generating 
functional for the Thirring model at finite T and fx = has been obtained in terms of the fermionic one (with an 
external electromagnetic source) for the Schwinger model . The Thirring model at nonvanishing T and fx has been 
analysed in ]ia] for real time and in H] in the torus. 

In this paper, we shall study, first, the Schwinger model in a medium at thcrmodynamical equilibrium, by introducing 
both the temperature and the fermion chemical potential fx. By considering a nonzero fx, we are able to study the 
system when there is a finite net fermion charge density (in the free case, fx is just the Fermi energy). Some of the 
questions that naturally arise are: i) whether one can provide simple solutions for the Schwinger model with a non zero 
fx, which extend previous studies non trivially, ii) whether bosonization takes place at finite fermion charge density 
and, if so, which is the boson mass, iii) which is the net fermionic charge of the resulting system, that is, whether 
the fermions are still confined into neutral mesons and iv) how does the chemical potential affect charge screening. 
We shall try to give answers to all of these questions. The second aim of this work is to provide an exact solution 
for the Thirring model at T ^ and fx ^ in the imaginary time formalism, in terms of the corresponding one for 
the Schwinger model, to compare with previous findings by other authors, using different methods, as a search for 
consistency and to get some new results. 

The plan of this paper is as follows. In section [H], we shall deal with the generating functional Z of the Schwinger 
model at nonzero T and fx, analysing several important items: the fermionic generating functional Zf with an external 
electromagnetic field, the role of the zero modes, the determinant of the Dirac operator, etc, by following steps similar 
to those in 0. The lack of hcrmiticity of the Dirac operator and a non-trivial phase factor J will be genuine and 
crucial features of the fi ^ case. They both will make necessary to extend the methods developed in H. From 
section III onwards, we shall restrict ourselves to the trivial sector, which is the only relevant one in order to study 
the thermodynamics of the system. We shall get Z, Zf and J by using functional methods, generalising what was 
done for fx = in B and deriving the proper extension of the point-splitting regularization when ft ^ 0. Section IV 
is devoted to: a) several physical results for the Schwinger model: the fermion charge density, the thermodynamical 
partition function, the boson propagator in the trivial sector, the Polyakov loop (the order parameter of the confining 
symmetry) and the screening length, b) the consistency of our methods. The tasks of obtaining an explicit solution 
and new results for the Thirring model at nonzero T and fx are undertaken in section Section ^ contains the 
conclusions and some discussions. Several results pertaining to the zero modes in the Schwinger model at nonzero T 
and fx are collected in an Appendix. 
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II. THE GENERATING FUNCTIONAL AT FINITE TEMPERATURE AND DENSITY 



Our starting point will be the generating functional for the Schwinger model in the imaginary time formalism of 
Thermal Field Theory [l5|Jl^| . We shall work in Euclidean two-dimensional space- (imaginary) time. In principle, we 
shall keep the length of the system L as finite, by imposing suitable boundary conditions in the spatial direction (see 
below). Thus, one properly defines the spectrum of the Dirac operator and avoids infrared divergences |7|Js||L4|] . At 
the end of the calculations we shall take the L — > oo limit. The finite density effects will be implemented by including 
a chemical potential [i associated to the conservation of the total electric charge (or the number of electrons minus 
that of positrons). Let A = (A 11 ) — (A , A 1 ) be the electromagnetic potential. Then, the generating functional reads 



Z[J,£ft=N(p) / VA exp 

J periodic 



d 2 x(T[A]+JA) 



T[A] 



1 



E 



1 

2^" 



where the fermionic generating functional is 



Zf[A,£,Z] 



(1) 



J antiper iodic 



d 2 x (-V P{A; n)il> + £,ip + V>£) 



and the Dirac operator is given by 



P(A;n) =@-ie4- ^7° 



(2) 



(3) 



In the above equations, N([3) is a temperature dependent normalisation constant, f3 = l/T, T being the temperature, 
J T is the integral over the Euclidean two-dimensional torus [0, 0\ x [0, L) and e is the electric charge, which has 
dimensions of energy. The fermionic and bosonic external sources are £, £ and J respectively. The electric field is 
E = Fqx = d{)Ai — 8\Aq and a is the covariant gauge-fixing parameter. It is important to remark here that the above 
covariant gauge fixing does not fix the gauge completely on the torus. There is still some residual gauge arbitrariness 
related to global gauge transformations, which we shall deal with later. The Faddeev-Popov determinant has been 
absorbed in the measure in (|l|), as it plays no dynamical role. Our conventions for the Euclidean Dirac matrices 
({7^,7i/} = are: 7 = 70, 7 1 = 71 and 7 s = — i7°7 1 are the Pauli matrices. 

The electromagnetic field and the bosonic external source are periodic in the Euclidean time with period /3 whereas 
the fermionic fields and sources are anti-periodic. An alternative approach, which we shall not follow here, would have 
been to take p(A;Q) in (||), with fermions satisfying the boundary condition tJj(x° + (3, x 1 ) = — e-Kp(Pfi)ip(x° , x 1 ). 
Concerning the spatial boundary conditions, they cannot be chosen as periodic, in general (after the above choice for 
the temporal ones), as the Dirac operator may have zero modes on the torus (to avoid duplications, we refer to R] for 
a justification). Without loss of generality, we shall choose A^ so that A M (x°, x 1 + L) — A^x , x 1 ) = d^—Qx^/eP) 
and hence 



+ L) = ip{x°, x 1 ) exp (i^ x ° 



(4) 



with $ the total flux of the electric field over the torus 



$ = e 



d 2 xE(x) = 2tt(ti + — ri_) 



(5) 



where n± are the number of zero modes with positive and negative chirality. The relation (||) follows directly from 
the axial anomaly 0,|l7]]. We shall introduce k = n + + ri_, the total number of zero modes. The gauge sector with 
k = will be referred to as the trivial sector. For later use, we recall the following factorisation property of the Dirac 
operator 



P(A; (j,) = exp(x°fi) ft (A; 0) exp(-a; ^) 



(6) 
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A. General structure of Zf[A,£,£] with zero modes 



The contribution of the zero modes to the generating functional has to be analysed carefully, in order to properly 
define the functional determinant of the Dirac operator. For that purpose, we shall follow the same steps as in M. 
However, there is an important distinctive feature of the fi ^ case, namely, that ip[A; fi] is non-hermitian. Hence, 
the set of eigenfunctions of ip is no longer an orthonormal basis in which the spinor fields could be expanded. To 
avoid this difficulty we shall expand the spinors in the basis of the hermitian operators P'P and pp^ . This will allow 
to separate the zero mode contribution up to a phase factor. We shall discuss below this factor and its relevance to 
the calculation. First, let us consider the set of eigenfunctions of the hermitian operators 



H(A-ix)<t>n= PHA;fi) p(A;fi) cj> n = fi n (f>, 



H(A; fi)tp n = p(A; fi) p\A; fi) ip n = ft n (p. 



(7) 



The operators H and H have the same eigenvalues /i„ > (for fi n > 0, p<j) n is an eigenstate of H) and the zero 
modes of H (H) are the same as those of P (P^)- In addition, since the anomaly (|j) is ^-independent (for general 
results on the independence of ano malie s on thermal effects, see n + — n_ is the same for both P(A;fi) and 

P(A; 0). As we shall see in section II C , all zero modes have always the same chirality. Therefore, the number k of 
zero modes is the same for H, H and P(A; 0). 

At this point let us expand the spinor fields ip and ip as 



p=l q=k+l 
k oo 



P =i 



q=k+l 



with 



q = k + 1, . . . , oo 



(8) 



(9) 



4>p (<^p) being the zero modes of H (H). In this basis, we have (ip q ,tp q >) — (<p q ,(/) q i) — 5 qq >, where (x,"0) = J T d 2 xx^ip 
is the scalar product on the torus. We get for the fermionic action 



d 2 xip P(A; fj,)ijj 



q=k+l 



[J<qf3 q f3 q 



(10) 



Then, the action is diagonal in this basis and, by doing the integration over the grassmanian variables a p , the 
contribution of the zero modes can be factorized. A crucial point should be noticed here. As the spinors tjj and ip are 
expanded in different basis, the Jacobian of the change of basis from T>i\) to Y[ p q da p dp q is not the inverse of that from 

T>tp to Yip q da p d(3 q . This fact was already noted by Fujikawa []l9f in the context of anomalies with non-hermitian 
Dirac operators. Since both changes of variables are formally unitary, when doing them simultaneously we are left 
with some phase factor exp[iJ'(A; fi)]. Notice that J {A] 0) = since then H — H = —p 2 (A; 0). Also, in principle, the 
phase factor is different for every k sector, a feature to be reminded by means of a superscript (k). Thus, performing 
the Gaussian Grassman integrals over dadad/3d/3 we get 



Z / [A,Z,t]=e X p[iJ W (A;n)]exp 



-i / d 2 xd 2 y£,(x)G(x,y,eA; fi)£(y) 
Jt 



n / d 2 xd 2 yax)MxHm(y)ydet'H(A-,fi) 



P =1 J T 



(11) 



where det' is the functional determinant when the zero modes are omitted (or factored out) and G(x, y, eA) = 
S^!Lfe+i ~^i= ( t ) q( x ) l p\(y) is the exact fermionic two-point function, satisfying the differential equation 



P(A; fi)G(x, y, eA; fi) = S (2 \x -y)-J2 M^tiv) 

P =i 



(12) 
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The second term on the right hand side of the above equation is the projector onto the H zero mode subspace. 
For simplicity, we have omitted a superscript (k) in both G(x, y, eA; p) and det'H(A; p). From ( |TT| ) we see that the 
zero mode contribution can be factorized in this basis, in which we obtain the contribution of | dct p\ = (det'i?) 1 / 2 . 
However, we have still to clarify which is the role of the phase factor exp(ij'): this will carried out in the next 
sections. Let us now remind how to obtain different quantities of physical interest from (pi]). If we are interested in 
the thermodynamics of the Schwinger model, the relevant quantity is the partition function Z = Z(0,0, 0), so that, 
from (fTy), only the trivial sector contributes: 



2(0,0,0) = N 



VAexp 

periodic,® 



iJ {Q) (A;p) + / d 2 xT[A] 



^detH(A; p) 



(13) 



We can obtain thermodynamic observables such as the free energy and the particle charge density, by differentiating Z 
with respect to the temperature and the chemical potential respectively, thereby generalising for finite charge density 
the study carried out in ||. The subscript '0' in the functional integral above indicates that only the trivial sector 
contributes. Then, in the trivial sector J can be identified with the phase of the fermionic determinant. We can 
also calculate the average fermion charge density p = L^ 1 J dx 1 (ipj°ip) in terms of the two point Green function. 
A remark is in order here: the equations of motion for the A field imply d\E(x) — —ieip'j ip, which is Gauss law. If 
the latter is imposed as a quantum constraint on physical states p0| , then it would imply p — with the boundary 
conditions chosen. As it is customary one may consider that an external compensating (say, ion) charge p ex is 
present to ensure charge neutrality and hence Gauss law holds for the total charge p to t — P + Pex = 0. Alternatively, 
one may consider an open system that exchanges particles with a reservoir ensuring charge neutrality. With this 
in mind, we shall make no further reference to p ex and concentrate only in the fermion charge density p for the 
electron-photon system. Therefore, (O) yields: 



iJ^(A;p) + [ d 2 xT[A 

JT 



i 1 
ZL 









/ VAexp 




J periodic,® 



x tr[7°G(a;, x)] +i 



VAexp 

periodic,! 



\/det H(A; p) 

iJ {1] {A;p) + J d 2 xT[A] v /det'iJ(A; A i)^ t 1 (x)7>i(a;) 



(14) 



and the k > 1 contributions vanish. We shall show in section II C that the zero modes of H and H have all the same 
chirality, given by sgn(<I>). Hence, <p\ and 4>i are both eigenstates of 7 5 with the same eigenvalue and therefore the 
second piece in the above equation vanishes. Then, only the trivial sector contributes to the fermion number density. 

We also remark that the property {j$,G(x,y)} = 0, which is not difficult to prove with the above definitions, 
implies, as in the case of finite T but vanishing p [Q, that the chiral condensate (ipP±i>), with P± — (1 ± 7 5 )/2, 
does not depend on G(x,y). However, from ( pd| ) we see that it will contain the phase factor exp^J^ 1 ) (A; p)] (see 
Appendix) . 



B. The imaginary part of the effective action 

As ip 7^ (ipy , we have found the extra factor J(A\ p), which is the source-independent piece of the phase of the 
generating functional. We shall analyse here its physical interpretation, at least in the trivial sector. The general 
form of J^{A; p) can be inferred from the symmetry transformation properties of the phase of the different quantities 
obtained from Zf after switching off the external sources. For instance, in the trivial sector, the object of interest is 
the effective action Zf[A, 0,0]. Now, recall that the /i-dependent term in the Dirac action is odd under the operation 
of charge conjugation C, since it is the number of particles minus the number of antiparticles operator. The rest of the 
Dirac action is even under C, so that C acts on Zf by replacing p — > —p' or, in other words, Z f [A c , 0,0] = Z) [A, 0,0] 
and, therefore, the phase of the effective action is odd under C, while the modulus is even. This is analogous to the 
case of the QCD effective chiral lagrangian, when the symmetry under consideration is spatial parity (P), the phase 
of the effective action being, then, the Wess-Zumino-Witten term p2fl . In our case, J^(A;p) should contain only 
C-odd combinations of the gauge field. As P is a symmetry of the effective action, J should be CP odd. In addition, 
it is not difficult to check that the p term does not generate any anomaly in the gauge current, so that imposing local 
gauge invariance (see comments below), the only term which fulfils such symmetry requirements is of the form: 
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J<® (A; fi) = F(T, fx,L) d 2 xA Q (x) (15) 



T 



where F(T, fi, L) is a function, undetermined so far (to be found explicitly later) such that F(T, fj,, L) = —F(T, — /i, L) 
since, by changing simultaneously /i — > — fi and particle by antiparticle, the theory remains unchanged. There is 
another point that is worth noticing here. Recall that in the torus the gauge transformations g (x°, x 1 ) : S 1 x S 1 — > U(l) 
are parametrised by TL x 2, corresponding to the two winding numbers (n, m) around the two circles. For any the 
most general gauge transformation A^ — ► A^ + d^a, which keeps d^A^ fixed (so that d^d^a = 0) and leaves unchanged 
the boundary conditions (in space and time) for both fermion and gauge fields is: a(a;o,a;i) = (2irnx )/ f3+(2Trmx 1 )/L^ 
up to an additive constant. Different choices of (n, m) correspond to nontrivial, homotopically disconnected, gauge 



transformations. But then, we note that the integral in (15) is precisely equal to n when A is a pure gauge field 



Hence, ( |15| ) changes by nF when we perform a gauge transformation g labelled by (n, m) and then it is not gauge 
invariant under nontrivial, gauge transformations. In this sense it is a topological term. Therefore, we are imposing 
local gauge invariance but still allowing a noninvariant topological term dependent on the chemical potential. This 
assumption is partially motivated by previous works in which /i-dependent topological effective actions were obtained 
p3LE4|. In the next section, we shall get F explicitly, through another method, thereby justifying our assumption. 



C. Instanton decomposition 

Following |t]], we shall decompose the gauge field into an instanton part A^ and a local fluctuation tfi as 

A/j, = An - e^dvcj) (16) 

with 4>(x) periodic in both space-time directions. Note that A^ yields a constant electric field E = $/(eL/?) and, 
hence, E = E + A<fr. By following steps similar to those in M, we separate first the contribution of A and </> in 



V det' H . Using that ft(A; 0) = exp (e7 5 0) ft (A; 0) exp (ej 5 (j>) and (§|), it follows immediately that 

p{A- n) - exp(e 7 5 0) p(A; fi) exp(e 7 5 0) (17) 

Notice that the operator H(A; fi) in (0) can be cast as 

H(A ]f i) = -V„V,(A,/i)- e ^ 7 5 

V„ = d u - ieA v + i^ b 8 ul (18) 

and the operator H is obtained from H by changing /i — > — fi. Hence, H(A;fi) — -V„V„(i;;i) — 7 5 -^, so that, as 

-V„V„ in (n8h is a positive operator, all the zero modes of H (A; fi) have the same chirality, equal to the sign of $ 
(recall that fy°, H] =0). On the other hand, from (|l7|), we get a zero mode of fl(A; fi) by multiplying a zero mode of 
P(A: ix) by cxp(— cy 5 <f)) = cxp(— esgn($)<f>), which is in turn a zero mode of H(A; /i). We can apply exactly the same 
argument to H. Then, the z ero m odes of H and H in (0) have both the same chirality, which is equal to sgn(<E>). This 



was already used in section II A , in order to omit the one zero mode contribution to the particle density. 

It is possible to separate the contribution of det' H (A; fj,) in det' H(A; /i), for arbitrary k. We have sketched the 
derivation in the Appendix, the general formula, for any k, being given in ( |A2| ). From that expression, we read the 
usual induced mass term for the boson field, with mass m = e/y^, which is independent on both the temperature 
and the chemical potential, thereby generalising the result for \i = previously derived in j^,f^| . Let us quote here the 
result for the trivial sector k — 0: 

e 2 f 

dctH(A; (j,) = dctH(A; /u.) exp — / d 2 x(t)(x)A(j)(x) (19) 

ti" Jt 



D. The determinant of the instanton operator 

In order to complete the analysis in the previous section, one should still study the spectrum of H (A; /i), which will 
be the purpose of the present section. First, by following 0, we shall decompose the field A as follows: 



G 



7 2tT, 



eL(3 (3 
ft A 



(20) 



which is the Hodge decomposition of the gauge field in the torus. The contributions proportional to ho and hi are 
the so called harmonic parts and are essential to correctly qu antis e the model [|7|Jl4|] . Notice that under a nontrivial 
gauge transformation (n,m) of the type commented in section [IB, the h fields above are the only ones changing and 
they do so as ho — » ho + n and h\ — *■ h\ + to, even for 4> = 0. The A-dependent terms in the last two equations are 
pure gauge contributions with A periodic in x° and x\, which will not play any physical role. For instance, with the 
covariant choice d^A 11 — 0, A is just a constant and that term does not appear in (p0|). Let us consider first the case 
$ = 0, which is the only relevant one for the partition function. 

Since 7 5 commutes with H , we choose the eigenfunctions of H(A; fi) as states of definite chirality, that is: 











(21) 



with 7 5 v]/ ± = ± 1 $ ± 1 since 7 s = diag(l, -1). Then, for $ = 0wc have to solve H ± (A; fj)^ = X ± (j) ± with 

ff ± (i; fi) = -(d - ih Q ) 2 - (ft - ihi ± ifi) 2 (22) 

and (anti-) periodic boundary conditions in the (time) space direction. In the above equation we have introduced 
ho = 27teho/[3 and hi = 2irehi/L. The eigenfunctions are plane waves and the corresponding eigenvalues are: 



1 (3 



n + - - eh 



27r 

——{k — eh\) ± fi 



(23) 



with n, k integers. Notice that the chemical potential breaks the chiral degeneracy which was originally present in the 
fi = case. Now, by using logdeti? = Trlogi/ and (23), we get: 



00 

logdetif(i;M) = r ^\og[(oj n ) 2 + {uj k -h 1 ± fi±iho) 2 } 



(24) 



±± 



with cj„ = (2n + l)n//3 and uik — 2-xk/L. Now, let us add and subtract log/3 2 J2 n k to * ne right-hand-side of the 
above expression. This procedure will give rise to a T-dependent infinite constant, which, in turn, will be absorbed, 
as customarily, in the normalisation constant N((3) [| 15 1^] . We can perform the summation over n in the above 
equation with the help of the two formulae p9,E5(| : 



log[(2n + l)V+/3 2 K±a) 2 ] 

OO 

E 



do 



i ¥+(2n + l) 2 ir 2 



1 



(2n+ l) 2 + 1 



log [l + (2n+l) 2 7r 2 ] 

i(i__jL_ 

6 \ 2 e 9 + l 



(25) 



In so doing, we obtain, finally: 



logdetH(A;fi) = ^ <j 2/3 (u> k - hi) + ^ log [l + exp -f3 (uj k - hi ± fi ± ho) 

±± 



(26) 



k— — c 



up to an irrelevant T and fi independent constant. In order to obtain the full partition function, we have to multiply the 
above expression (which for e = reproduces the partition function for free fermions at finite density) by exp(ij^) 
in ( |13j ) and by the ^-dependent one in ( |l9| ) and, then, integrate over the gauge fields (</>, ho, hi). As a consequence of 
the decomposition of the gauge field chosen here, the phase factor only depends on ho, so that: 



f d 2 xA (x) = h 2irh L 

Jt 2e 



(27) 
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since (j> is periodic in the space direction. Then, the <j> contribution to the partition function in (|19|) gives the partition 
function of a free massive boson 0. All the dependence on the chemical potential is included in the (ho, hi) part, as 



given in ( |26| ) and (27). However, we have still to determine the value of F(T, H,L) in (jlq). Before undertaking that 
task, and for completeness, let us analyse the spectrum of H(A; /x) when $ ^ 0. In this case, we have 

^(A; n) = -{do - iho + ^xif ~ (ft ~ %h~x ± *m) 2 T eE (28) 

together with the boundary conditions in the spatial direction given in (|]). This eigenvalue problem is solved in the 
Appendix. From the result foun d th ere, we remark here that the fi dependence, when <!> ^ 0, appears only in the 
s tat es, while the determinant in (A7) depend s on the temperature T but not on \i. As the norm of the zero modes in 
(A8) is also fj,- independent, then, if we go to (|a|) and (|ll|), we realize that the dependence of the fermionic generating 
functional on /i, when $ ^ 0, is encoded in the determinants of the matrices in ( |A2j ), that follow immediately from 
the spectrum found for a = 0. Besides, there are fi dependences in both G(x,y, eA;fi) and (A; fx) . 



III. THE GENERATING FUNCTIONAL IN THE TRIVIAL SECTOR 

In this section, we shall use functional methods in order to calculate the generating functional in the trivial sector. 
By employing these methods, we shall also obtain the phase factor J^\A\ /i). This will allow to get the full fermion 
charge density and the partition function, as well as to establish the consistency with the results of the previous 
section. We remind that the covariant gauge-fixing is not complete on the torus. We still have the freedom of 
performing a nontrivial gauge transformation (ho, hi) — * (ho + n, hi + m), with the h fields in ( |2C| ) and n, m integers, 
corresponding to loops that wind n times around the temporal direction and m times around the spatial one. It is 
clear that n and m are not fixed by d^A^ = 0. Throughout this section, though, we shall work with the covariant 
gauge-fixing, ignoring this residual gauge arbitrariness. The latter has also been treated in a situation related to 
the one analysed here, but not quite identical with it: specifically, for (real time) QED on a spatial circle, at zero 
temperature and chemical potential p6fl . We have to bear in mind that in ([!]), we are integrating the gauge field 
over all possible values of the fields h^, that is, /i M € IR. Fixing the gauge for those fields would consist in restricting 
them to a range [0, 1] x 7L since they change by an integer under a global gauge transformation. Then, if the 
effective action is globally gauge invariant, and unaffected by the residual gauge arbitrariness, the difference between 
integrating over all h or restricting them to a [0, 1] interval, is an infinite constant independent on T and /i. So that 
arbitrariness cannot affect physical observables such as the free energy or the particle density. For e = 0, the action 
depends on derivatives of the electromagnetic field, and the covariant gauge fixing, even if not complete, suffices to 
get a well-defined propagator, unaffected by that arbitrariness. In general, when e is non- vanishing, Z[J, is also 
unaffected by the arbitrariness, after having integrated over all fields. However, for a given A„, both det H(A; /i) and 
the fermio nic g enerating functional Zf may be subject to it, even if [i = 0. In particular, when [i ^ 0, we have seen 
in section |HBj that Z$ is not globally gauge invariant, due to the induced topological term in the phase J, which 
changes when h — > h + n. Thus, if we restrict h to a [0, 1] interval, the result for the observables woul d dep end on 
our choice and then it is consistent to let ho £ JR. We shall come again to this point at the end of section [V A , where 
we shall perform the integration over the h fields explicitly, using the results derived in section [□], It is not difficult 
to check that all the formal functional manipulations that we shall carry out in this section, except those related to 
L[A], are also unaffected by the residual gauge arbitrariness. 

Thus, as a first step, let us rewrite the generating functional in (Q) for the trivial sector, with the aid of standard 
functional techniques, as: 



Z [J, £] = Z E mZ f exp 



x exp / d 2 xd 2 y 

IT 



le I d x - , , 7^- — — 
l -J^{x)D^(x - y)J v {y) - i£(x)S(x, y; ^(y) 



(29) 



where Zem and Zf (Zf(T, /i, L) = Zf for short) are the free boson and fermion partition functions and D^ix — y) 
and S(x,y;fj.) are the free gauge boson and fermion propagators respectively: 



D llv {x-y) = — , 

7i,k=—oo 



1 



X i ( ^s^^Wv 
<V + \ a ~ l ) 9~ 



S{x,y;n) 



(3L 



E 



n,K— — oo 



7°(w„ + ifi) + 



(30) 
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where lu — (cj n , u>k), with u>k — 2irk/L, lu h = 2nn/P in the bosonic propagator and w„ = (2n + l)it//3 in the fermionic 
one. Now, we shall make use of some known functional differentiation formulae 113] , and, in particular of: 



exp 



,2 , 2 8 .5 
a xd y —. . A(x,y) 



= exp 



xexpi J d 2 xd 2 y£(x)B{x,y)£{y) 
d 2 xd 2 yl{x)B{x,y)£,{y) + L 



(31) 



Here, A(x, y) and B(x, y) are arbitrary functions, to be regarded as the kernels of the operators A and B, respectively, 
B = B(l + AB)~ X and L = — Trlog[l + AB}^ 1 , Tr indicating the trace over functional and Dirac spaces. Thus, one 
finds: 



Z[J,£,l\=Z EM Z F exp 
1 

x exp -- 



f 
2 

d 2 xd 2 y 



d 2 xd 2 yJ tl (x)D^{x-y)JM 



D^(x-y) 



SA^x) 



x exp <^ -i / d 2 xd 2 y£,(x)G{x 1 y,ieA;n)£ > (y) + L[A] 



with An(x) 



(32) 



; J T d 'ryD^^x — y)J v (y), after having performed the functional differentiations, which appears to 



leave no trace of the residual gauge arbitrariness in Z[J,£,Q. The so-called closed fermion loop functional L[A] can 
be written formally as: 



L[A]=ti D / de' I d 2 x A'(x)G(x,x,ie'A;[i) 



(33) 



where tr^ denotes the Dirac trace. We remind that G(x, y, eA; fj,) is the two-point function, which, in the trivial 
sector, fulfils Eq. ([l2|) with k = 0, that is, with its second term on the right hand side omitted. 

In order to get a well defined expression for the generating functional, we need, first, to regularise L[A] in (^3|). For 
that purpose, we shall appeal to the point-splitting regularization fl2]| . Therefore, we should deal with the following 
limit: lim^^j, G(x, y, eA; fi). In Minkowski space-time, the limit should be taken by keeping the points x and y 
relatively space-like, in order to maintain causality jl2|. As we are working in Euclidean space-time, we shall not 
impose, in principle, such a restriction. We shall comment below on the different ways of taking the limit. Before that, 
and generalising |fl2| , we shall derive the point-splitting regularization prescription in our present case with nonzero 
fi. We start with the formal definition of the gauge current in the presence of an external background field A^: 



{3^{x))f[eA] = (^(x)7^(a;)>/[eA] =i lim tr D 7 M G(a;, y, eA; fi) 



(34) 



where (O) f = j VipUipO exp[— J ippip}. We shall obtain the regularised version of the right-hand-side of the above 
equation as follows: we shall demand that such a regularised gauge current be conserved and gauge invariant. Notice 
that, under a gauge transformation A^ — > A^ — d^K, G(x, y, eA; /i) changes as: 



G(x, y, eA; fi) -> G(x, y, eA; n) exp ie[A(ai) - A(y)} 
Based upon this, it is easy to show that the product: 



G(x,?;,eA; M )expHe P ' d£° A^)] 

J X 



(35) 



(36) 



is gauge invariant. However, if we replace (36) into (p4[), calculate the divergence of the current so defined, and 
use p(A; fi)G = S^ 2 \x — y), we find that such a divergence does not vanish for fj, ^ 0. To ensure that the current is 
divergenceless, we have to add an extra /i-dependent term, which leads to the regularised gauge current: 



{^{x)) r f e9 [eA] =ilimti D j IM G{x,y > eA;n)exp[-ie / d£ a 4,(0] exp[-/x(a; - y )] 

J x — m I 



(37) 



which is, indeed, gauge invariant and divergenceless. Note that Euclidean covariance is broken since the system is in 
a thermal bath. We are now ready to define the regularised fermion closed loop as: 
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L r ^[A] = -itr D J* de' J^xA^j^Y^lie'A] 

= tin / de' / d 2 xAn(x) lim tvolu.G{x, y, ie'A; n) 

Jo Jt x ^y 

x exp[e' f d^MO\ exp[-/i(x° - V°)} 



(38) 



The limit x — ► y has to be taken in a symmetric way, regarding (x, y) ]T^ |. In order to calculate the fermion closed 
loop in (|3^), we shall consider an ansatz for the exact Green function similar to that in ||: 



G(x, y, eA; /i) = exp[— ie\x(x) - x(y)]S(x, y; /i) 



(39) 



It is not difficult to check that, with the above ansatz, G(x,y) is a solution of pG(x,y) = 8^ 2 \x ~ y) , provided 
that x( x ) be a solution of ft x x(x) = —^L(x). In turn, the solution is: 



X(x) 
A(x - y) 



d 2 yA(x - y) @ y /(y) 



^ CO 

w. E 



•(»-!/). 



/3L 



(40) 



,k— — c 



where a; = (w„,Wfc) with w„ = 27m//3, = 2nk/L. So G(x,y, eA; fi) appears to be unaffected by the residual gauge 
arbitrariness. 

Hence, from the above equations, we have to find the behaviour of S(x,y;/i) when x — > y. For that purpose, we 
apply the standard formulae: 



oo 



•i/z) 



1 

2^ 



dwf{uj) + d> duofiuj) 
Jc 



oo=p?e+i/i 

du)f(u) 



dujf{uj) 



1 



j±i/8w _|_ ^ 



dujf(u)- 



1 



oo=p*e 



(41) 



where e — ► + , w„ = (2n + l)m//3, LOk = 2nk/L and C is the rectangular contour in the complex plane running 
through the points (+oo, — oo, — oo + i/i, +oo + z/i). We apply ( |4l|) to ^(a;, y; /i) in ( [j0| ) and retain only the dominant 
contributions in x — > y. The complex plane integrals along the lines (— oo=pze, +oo=Fie) and (— oc^fie+ifi, +oo^fie+i^i) 
are performed by forming a closed contour, by means of an infinite arc above (below), corresponding to the +(-) sign 
in (41), and applying the Residue Theorem. We get: 



lim S(x,y;n) 



2tt [x — yY 



where 



1 f°°. 

F ^,L) = - £ 



k— — oo 



l 



l 



(42) 



(43) 



with LOk — litk/L. We recall that (|42j)reproduces the T 
{Fiji = 0) = 0) in § 



e /3(^fc+M) _|_ 1 e /3(wfc-/j) _|_ I 

H = result given in Il2| and the /i = 0, T / one 



Next, we shall replace both ( |f0| ) and the limit (42) into (39) and (p8[). We have taken the x — ► y limit in two 
different ways and established that the same result is arrived at. We have taken, firstly, x° — y° — ► 0, x 1 — y 1 — ► with 
(x 1 — y 1 )/(x° — y°) = 1 and, secondly, the Minkowski causal choice (see |l2|) x° — y° and (x 1 — y 1 ) — > 0. Anyway, 
what it is important to note here is that the exponential \x dependence in (|42|) is exactly cancelled with that in the 
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regulator in (|38|). Then, no matter how we take the x — > y limit, we always get a term tr J ^(x)^f°F in L reg . The 
possible divergence in L reg arising from the first piece in (|42j) is absent since we have taken the limit symmetrically. 
Thus, finally we arrive at 



1 



L reg [A} = - / d 2 xd 2 yA a (x)U a ^x-y)A^y) + 2eF(T,^L) / d 2 xA (x) 



with 



1 e 2 



Sal 



UJ a UJf3 



(44) 



(45) 



where (jj„ = 2im/(3 and = 2nk/L. So, L reg [A] is, in principle, affected by the residual gauge arbitrariness. The 
functional differentiations with respect to A(a;) in ( p2\ j can be performed by employing the following formula, which 
is valid for any linear operators P and Q fl2l 



(46) 















exp 


■_if± p ±- 
2JSA 5A_ 


exp 


l -jAQA + iJf-A 


= exp 





if A(l-QP)- 1 -f+±tr\og(l-QP)- 1 + l - J /Q(l-QP)- 1 / 



where P = P(l — QP)^ 1 and we have omitted, for simplicity, all the space-time dependences. We shall concentrate 
on the bosonic generating functional Z[J, 0, 0], as we have already analysed the dependence on the fermionic sources 
in the previous section, up to the phase factor (to be calculated below). Upon applying ( |4^ ) in (|3^), L[A] being given 
in (44), we get: 



Z[ J, 0, 0] = ZemZ f exp 
1 



x exp 



d 2 xd 2 y(J - iG),(x)Vf u {x - y)(J - iG) v (y) 
Trlog(l + ILD) -1 



(47) 



where Go = 2eF(T, /x, L), G\ — 0. In turn, ID^] is the exact boson propagator at fi = 0, which can be expressed, 
formally, as: E> (0) = D(l + WD)' 1 , D and n being given in d3Q) and (Eq), respectively. Its explicit representation is: 



oo 



n.k— — oo 



i 



Sfiv - 



(48) 



where m 2 = e 2 /ir is the induced boson mass. 

In the following section, we shall obtain, from ( [i"4| ) and (|47]), the complete form of the fermionic generating functional 
in the trivial sector, including the phase factor, the exact boson propagator at fj, ^ and the partition function and 
check the results with the method used in the previous section. 



IV. PHYSICAL RESULTS 



A. The charge density and the partition function 

By setting J = in the expression ( p!7| ) we get the partition function 

Z(T, /i) = 7 y ,fl) / (T, /i = 0) exp - \ [2eF(T, n,L)] 2 f tfxfyto® (x - y) (49) 

However, the above integrals of Eoo are ambiguous, in the sense that changes in the order in which such integrals 
are done yield different results. For instance, let us perform, firstly, the spatial integrals: that would force us to set 
k = in (|48|), which would give rise, after doing the temporal integral, to an infinite and gauge dependent result. 
Conversely, by changing the order of the integrals, let us perform, first, the one over the (imaginary) time. In so 
doing, we arrive at a finite and gauge-independent answer. The latter prescription seems to be a natural and physically 
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reasonable choice. In addition, as we shall see, it is consistent with the results obtained in the previous section. If we 
adopt this prescription in (p!9|), we get: 



Z(T,p) 



Zf{T, p) 
Z F (T,n = 0) 



Z(T, (1 = 0) cxp[-2/3LirF 2 (T, p, L)\ 



(50) 



This is a genuine nonperturbative result, since the argument of the exponential in the term that corrects the free 
fermionic partition function is independent of the electric charge. As we shall see, such a nonperturbative behaviour 
comes directly from the topological structure discussed in section namely, from the phase J\A\ p] of the fermionic 
generating functional. At this point, we recall the result obtained in for Z(T,p — 0). The latter partition function 
was proven to factorise into the product of that for a free fermionic field times that for a free massive boson field 
with mass m, divided by the one for a free massless field (finite temperature bosonization). In our present case, with 
p =/= 0, we may wonder if such a factorisation actually takes place as well and, if so, whether the whole system may 
have a net fermionic charge or not. Let us consider, first, the free fermionic partition function Z F (T,p): 



log Z F (T,p,L) = 



k= — o 



{log 



a -/9(wfc-/i) 



log 



I _|_ e -P{uk+tA 



The net free fermion charge density is 



±~ log Z F (T, p, L) = -2F(T, (i, L) 



F being given in (43). Thus, from (p2J) and (EG) we have 



Z{T, p) = Z{T, n = 0,L) cxp <^ -2j3L 



nF 2 {T 1 p,L) + [ dp'F(T, p', L) 
Jo 



(51) 



(52) 



(53) 



which is our final result for the partition function at /i ^ for finite L. Now, let us analyse the behaviour of the 
function F in the limit L — > oo. In such a limit, the sum over k becomes a trivial integral, which yields: 



F(T,p,L^ oo) = - — 



(54) 



Thus, in the L — > oo limit, the p dependence in (p3) exactly cancels out. We get for the full partition function: 



Z(T, p,L—> oo) = Z(T, p = 0, L — > oo) = ZemZ f (T, jj, = 0, L — * oo) 



x exp 



^Trlog(l + ILD) -1 



(55) 



which is only T dependent. Its explicit expression can be found in section IV of Q. That is, in the L — > oo limit the 
system bosonizes as well, and the only effective degrees of freedom are those of a massive boson field, the net fermionic 
charge of the system being p = ((3L)~ 1 d\ogZ / dp = 0. We remark that: i) This is a nonperturbative effect (which 
has been established due to the fact that the Schwinger model at finite chemical potential and temperature can, still, 
be solved exactly); ii) It is characteristic of two dimensions. Recall for instance that in pcrturbative four-dimensional 
QED in the infinite volume limit, the free energy \ogZ is /i-dependent and hence p ^ |pl| ; iii) The result fl55|) holds 
in the L — > oo limit: if we keep L finite, F is no longer given by ( p4[ ) but it acquires further corrections and the study 
of the counterpart of ([55]) will not be attempted here. 

We now turn to the issue of the C-violating phase factor, addressed previously in section ||. Thus, we consider 
the fermionic generating functional Zf\A . £, £1, treating now A as an external background field. We apply (31) and 
we arrive at an expression analogous to (p2[), but now taking J = 0, omitting the derivatives with respect to A and 
replacing e — » —ie. By recalling the regularised fermion closed loop functional obtained in (B80, we get: 



z f [A C> I] = Z f Z E m exp 



-i / d 2 xd 2 y£(x)G(x,y,eA;p)£(y) 
Jt 



x exp 
x exp 



d 2 xd 2 yA a (x)IL^(x-y)A (y) 



-2ieF(T, p, L) / d 2 xA (x) 
Jt 



(56) 
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The comparison with (|ll]) in the trivial sector leads to identify: 

jW(A,n) = -2eF(T,ii,L) J d 2 xA (x) 



(57) 



that has the form that we had anticipated in (|15|), based upon the C symmetry properties of the phase of the fermionic 
determinant, and, so, F = —2eF. 

We shall now provide an interesting check of consistency. On the basis of ([l3|), (|l9|), ( p6[ ) and (j27|), we can calculate 
the partition function through the method developed in section ||. We shall do so in the L — > oo limit. Let us first 
differentiate with respect to fi in ( |26| ) and then take the L — ► oo limit by replacing ujk by continuous u> and ^ fc by 
(L/2n) J du>. The resulting integral can be done using 



+ OC 



du> 



1 



1 



1 + £0(i>J — a+b) 2_ 4- g/9(o)— o— 6) 



= -26 



(58) 



which is convergent when both pieces of the integrand are added together. Then, in the L — *■ oo limit we obtain: 



— logdet/f = 



and, hence: 



Vdet /x, T) = v / ddTff(l; M = 0, T) exp ( L/3— 

27T 



(59) 



(60) 



We now insert the above result in (|T^) and, then, in ([Tl|), for the trivial sector, and we compare to (^6|). We see 
that, before integrating out the photon field, the results obtained with both methods are consistent with each other 
in L — ► oo, since we get the same /i dependence, namely, the one coming from the phase factor (the explicit form of 
which has been obtained with the second method) times that of the free fermionic partition function, as it appears 
in ([33|) when F — > — /j/27r. To complete the check, we shall see that the integration over the gauge boson fields also 
gives the same result with both methods in the L — > oo limit, namely, that in (|55|). From our previous discussion, 
it follows that the only part that depends on fj, in the full partition function is the integral over the (ho, hi) fields, 

the integrand of which is the exponential of ( |27j ) for $ = times (— 2ieF), multiplied by J det H(A). From (|60|), it 
follows that all the h dependence of det H is contained in the [i — part. In order to extract the dependence on /to, 
we differentiate again in (E9). Using again (pq) we get in the L — * oo limit: 



_d_ 



-logdet/f(l;/i = 0;/i ,/ii) =2ie V / 



duj 



a=±' 

%irLe 2 h 

i9 



(61) 



which turns out to be independent on h\. Thus, from ( |60| ) and (^Tj), we derive: 



\/detif(I;/i,r; ho, hi) = Vdet if (I; yt = 0,T;0,hi) 



x exp ( Lj3 



2tt 



2irLe 2 hl 



(62) 



At this stage, we have to integrate over the fields (ho, hi). As commented at the beginning of sectio n it is 
necessary to integrate over ho € IR, to achieve consistency. This is reinforced by the topological term (|57|) in the 
phase of the fermionic generating functional. Thus, from ( p7|) and (|62|), it follows that the relevant factor carrying 
the [i dependence is: 



dho exp —2i[ihoLi — 

-oo V P 



= e\ — exp 
V 2L V 



H 2 Lf3 
2tt 



(63) 



It turns out that the exponential fi dependence in (63) cancels that in (|62J), and, hence, we arrive, again at the result 
(l55|), obtained with our previous method. 
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B. The boson propagator, the screening length and the Polyakov loop 



The exact boson propagator, that results after the integration of both fermion and gauge fields, can be obtained, 
in the trivial sector, just by differentiating (^) with respect to the external sources and setting J = 0. We find: 



TD^(x -y) = Djx - y) - (2eF) 2 [ d 2 ud 2 zTD^(a 

JT 



u)TD$(y-z) 



(64) 



D^(x — y) being given in (|48|). The above expression is, again, formal, in the sense that we have to specify the order 
in which the spatial and temporal integ rals s hould be done in the second piece on the right-hand-side. If we adopt 
the same prescription as that in section IV A (that is, by doing the integrals in the order that gives a finite answer) 
we arrive in momentum space at: 



D^(o; n ,cjfe) = T>^(uj n ,Lu k ) - 5 n0 6 k0 (3L 



47r 4 F 2 (7>,L) 



(65) 

Dqq^ [u> n , uJk) being given by the expression between square brackets in ([f8|). In the L — > oo limit, u k is replaced 
by a continuous variable to and L5koF 2 goes to (/i 2 /47r 2 )(5(a;), in (|65|). As a check of consistency, we notice that the 
propagator ( |65| ) (which is /i dependent) satisfies the gauge invariance condition w' i w 1 'D (u/ = a, with a the gauge-fixing 
parameter. For finite L, the only Euclidean pole of the propagator ( |65| ) remains at uj 2 = —m 2 — —e 2 /n, which defines 
the mass of the boson and is independent on /i. However, in the L — > oo limit, a new pole appears at u 2 — due to 
the contribution of S(uj) in the second (p, dependent) term. To clarify what this pole means physically, let us calculate 
first the inverse correlation length squared M 2 in the L — ► oo limit, through the usual definition 



M 



P> oo (w„ = 0, w -> 0) 



where T> a p(uj ni oj) is the boson self-energy, which is defined in momentum space through 

Accordingly, from ( |66| ) and the propagators in ( p^p| ) and (^), we obtain 

Af 2 (T,^) = 



(66) 



(67) 



(68) 



That is, the mass M 2 vanishes in the L — > oo limit. This is indeed a consistent result if we recall the connection 
between the screening length and the equation of state of the system given in |16|| : 



M 2 (T,ri=e 2 —p(T,[i) 



(69) 



p being the total fermion charge density. Hence (|68| ) and ( p9|) are consistent in the L — >■ oo limit: we got a total zero 
charge density, while M 2 tends to zero. However, were M 2 be interpreted as a vanishing screening mass, then the 
system would be in a confined phase and the Z symmetry would be restored, which does not occur for p = and 
T =/= in the L — > oo limit JlO) ] . We shall outline below a correct understanding of the screening mass. The order 
parameter of this symmetry is the thermal average of the Polyakov loop 



exp 



ie I drAo(r, x 1 ) 
o 



(70) 



for an arbitrary charge e and a spatial point x . In addition, the correlator of two Polyakov loops with charges e and 
— e and spatial points x 1 and u 1 measures, at large spatial separation, the effective potential between both charges. 
Thus, in order to clarify the role of the two different masses m and M, let us calculate the Polyakov loop and the 
correlator for fi ^ in the L — > oo limit. First, we shall use the functional methods developed in section 1H. After 
the obvious relabelling of variables, the insertion of Pe{x l ) in the functional integral (pi) amounts formally to replace 
J(y) — * iP(y) w ith Po(y) — ^{y l " xl ) an( i Pi(y) = 0. Then, recalling our result (|^we get 
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exp 



d 2 yd 2 zlD$(y - z)P (y) (P (z) - 2G (z)) 



= exp 



0j? 
2m 



(71) 



where in the last step we have first followed our convention of performing the time integration before the spatial one 
and then we have taken the limit L — > oo, so as to transform spatial sums into integrals. The result ( |7l| ) can also 
be obtained using the decomposition ( |l6| ) and (^) for the gauge field. As we saw in previous sections, all the \x 
dependence is included in the ft, -dependent part. Then, since 



drAoi^x 1 ) = 2irh -dl / dr</>(r, x 1 ), 



(72) 



all we have to do is to insert a piece exp(27ri/i e/e) in the integrand in (|63|), in the L —* oo limit. It is straightforward 
to get again the \x dependence of the Polyakov loop in ([n]) consistently, once we integrate over ho G H. 

We see that (Pg(a; 1 )) never vanishes for any value of T and fi, so that the Z symmetry is never restored for L — > oo. 
Our result extends that obtained in j[o| for fj, = 0. Then, there should exist a nonvanishing screening mass, which, 
by consistency with our result for the propagator, should be m — e/i/n, independent on fj,. This is confirmed by 
calculating the correlator of two Polyakov loops. Following the same steps as for (P) we obtain with both methods 
that (Pe(x 1 )P-e(u 1 )) is independent on \i and 



lim (Pg(a; 1 )P_g(u 1 )) 



exp 



0& 
2m 



(73) 



that is, m is the screening mass between the two opposite charges. It seems clear that the definition (pq), leading to 
(|68|), does not give the right answer for the screening mass between two opposite charges, nor can be used to infer 
any conclusion about the breaking of the confinement or Z symmetry. However, it is the consistent equation to use 
if we want to get the equation of state of the system through (p9|). 

Notice that we have restricted ourselves to the trivial sector. The boson propagator could receive additional 
contributions from sectors with fc > 1, apart from that in (64). However, we have seen that p only receives contributions 
from the trivial sector, so that, (68) is valid for any sector if (39) holds. In addition, it is enough to restrict ourselves to 
the trivial sector to calculate correlators of Polyakov loops [10|. Hence, we expect the above results for the screening 
mass to remain valid for $ ^ 0. 

If the length of the system L is kept finite, the definition (|66|) for M 2 is no longer valid. Using (^39j) instead would 
give rise to a nonzero and ^t-dependent M 2 , directly from (|33|). Remember that for finite L there is no pole of the 
propagator at ui 2 — 0. However, it is not clear whether (|6^) remains valid for finite L. On the other hand, for finite 
L and = 0, the screening mass is only different from zero if e/e £ 7L. [ fTpf . A more rigorous analysis of the finite 
length corrections for /x ^ is beyond the scope of this work. 



V. THE THIRRING MODEL AT FINITE T AND /j, 
A. The generating functional 



We shall consider a system of many massless fermions (and antifcrmions) in one dimension (inside a finite interval 
of length L) at equilibrium at absolute temperature T and chemical potential /i. By assumption, their dynamics is 
described by the Thirring quartic lagrangian. Let £, £ be fermionic external sources. Then, in the imaginary time 
formalism, the generating functional of the system reads now 



Z[^=N(0,[i) PV^exp 

J antiperiodic 

\9 2 



d 2 x 



/i 7 )^ + & + ^) 



(74) 



where g is the coupling constant. The partition function is Z[Q, 0]. One can also cast as follows 

5 



, 2 _J_ _S 

t X '^{x) lv SJ v {x) df( a 



x Zi[£,£;J]U =0 



(75) 
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Zi[£,&J\=NZ F exp 
x exp 



T 

d 2 x"d 2 y"l{x")S(x",y"-^(y") 

T 

+ 00 



iaj-(x'-3/') 



n,/c=— 00 



^ - f{u> 2 ) 



(76) 
(77) 



with S in (pO[), Z F the free fermionic partition function and J = (Jo, J\) a boson source, to be set equal to zero after 
all functional differentiations with respect to it have been performed in ( |75| ) and u> in ( [77| ) is the same as in the gauge 
boson free propagator in (|30|). We have introduced an arbitrary function fjui 2 ), by virtue of the fact that the current 
-07 1/ -0 in the actual Thirring model is conserved. A proof of equations ( J7q ) and (^) follows readily through steps 
similar to those in (fjjl^]. At this stage, using standard techniques |l2], one can rewrite eqs.([75|) and (|76| ) as 



exp 



M^(x') 



x Z f [A:,^^]\ A , = o 

i(J 



Z f [A* u ,£,Z]=NZ F ex V 
x exp 



d 2 x-J^Al(x)^- 



(78) 

(79) 
(80) 

Again, one sets A* = in the above equations, after having carried out all functional differentiations. Standard 
functional techniques allow us now to establish that Z/L4*,£,£], as given in (|7|) also coincides with the right-hand- 
side of (||) (when due care is taken of the normalisation factor N((3, [i)) provided that, in the latter, one replaces e by 
g and A by A*. Let us concentrate on A* belonging to the trivial sector (see comments below). Then, by recalling 
the developments in section [II, one finds immediately 



A* u {x) = -i d 2 zK uP (x ~ z)jP{z) 



d 2 x 



d 2 yax)G(x,y,gA*;^ay)+L[A* 



(81) 



where G(x, y, gA*; /1) and LL4*] are now given by the right-hand-sides of (^9|), ( [io|) and ( fl4| ) (with the same S, A, IT 
and F), when one replaces e,A by g,A*, respectively. Thus, we have provided the solution for the Thirring model 
at finite T, fi in terms of the fermionic generating functional for the Schwinger model, without zero modes ($ = 0). 
The formal use of eqs. 

©"© and ©"© ( which _ 

can be checked upon comparing the corresponding expansions 
into powers of g) is justified if we restrict Z /L4*,£,£] to the trivial sector of the Schwinger model, as implemented 
through (|sij). In this regard, it is interesting to note that, upon using ( ]8C| ) in (||) and taking into account that the 
propagator K a p in (77) satisfies periodic boundary conditions in both arand ar directions, one immediately finds 
<I>[j4*] = 0, which confirms that Z f[A* . £, £] should be restricted to the trivial sector when use is made of (|7^)-(|80|) 
and, hence, the consistency of (|8l|). On the other hand, this appears also to be consistent with the idea that, in the 
end, we are going to set A* = and then we can take the vector field A* v as a configuration in the trivial sector, that 
is, topologically connected with A* = 0. It is unknown to us whether the Thirring model may have other solutions 
(besides that given in ([7q) and (81)). 



B. The thermodynamical partition function and fermion correlation function 

The thermodynamical partition function becomes, upon applying (|i^) to ([78]) and ( |si| ) 

Z[0,0] = Z F exp I^Tilog(l + UK)- 1 
- \{2gF{T, /1, L)f j d 2 xd 2 y [K(l + UK)- 1 ] °° (x - y) 



(82) 
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where in momentum space we have 



, +00 



n,k— — 00 
+00 



[K ( i + n^)-T /3 (^-2/) = ^ E 



•(3-1/) 



•(a:-!/) 



sra/3 



1 + <7 2 /tt 



n,k= — 00 



1 + .g 2 /^ 



.9 2 A 
1 + ,g 2 /7r 



(83) 



Using (B4j), (pgh yields readily for L — * 00 



Z[0, 0] = Zp(T, fi) exp <j La log ; X 9/ _ \ exp j -/3L- ff ' 



+ g 2 /n 



2tt 2 



1 + .g 2 /Tr 



b f(0) 



g 2 A 

1 + ff 2 /7T 



(84) 



where 



+00 



dk 
2^ 



2 ^ 

n= — c 

& = lim / d 2 xd 2 y 

^ > °° (/?£) 2 7t 



•(»-!/) i 



(85) 



n,k= — 00 



Notice that the first exponential on the right-hand-side of ( |84| ) is independent on /?, /1 and then it is irrelevant as 
far as the thermodynamics of the model is concerned. On the other hand, the constant b (which is independent on 
both T and [i) could in principle give rise to a dependence on /(0). Ag ain, the very definition of b is ambiguous. If we 
agree to evaluate the (imaginary time) integrals over a; and y° in (|85| ) before the spatial ones, then 6 = and, hence, 
the independence of the partition function on f(u> 2 ) follows, which is a welcomed result. Another independent reason 
to favour this prescription is that it is the same as that leading from (^) to (|50|). Now, by recalling the expression 
for the free charge density in ( |52| ) , we have from (J8J) that the total fermion number density of the system reads in 
the L — > 00 limit 



P = 



(86) 



Hence we obtain that the Thirring model at finite T and fi is no longer a free fermion gas, but the fermion density 
acquires a correction in g 2 , as it stands in (|S6"1). It differs from the result in ||l3f in which only the free contribution 
to the fermion density remains. It is clear from our analysis starting from the Schwinger model that the correction 
to the free gas comes entirely from the topological contribution depending on the F function. This contribution only 
depends on the harmonic field ho in the decomposition of t he g auge field. The calculation in pl| was done in real 
time formalism, in which this term is not present, whereas in [jl4| the model is solved in the torus. As it is emphasised 
in plplO], the toroidal compactification is very useful to deal with infrared divergences and the harmonic parts of 
the gauge field are essential to correctly quantise the model. It is the most natural choice when using the imaginary 
time formalism, as we have done in this work. On the other hand, if we evaluate the pressure of the system, which 
follows directly by taking the logarithm of the partition function in (pi|), our result (not quoted for brevity) agrees 
with JTil , which provides a check of consistency between our methods and those used in that work. 

Finally, we shall give the exact fermion correlation function for the Thirring model at nonzero T and /1 



<5 2 log 
G{x,y) = — =■ 



= Q(x,y)S{x,y) 



c=e=o 

+00 



n,k— — 00 



g 2 A 

1 + <? 2 /tt 



2gF{T,n,L)-c 



(87) 
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where we have used again ( |4q ) into (78) and ( pl| ) and performed the functional differentiation. In turn, c is given by 
the formal expression 



0,k 



+ /(w ) 



1 + g 2 /ir 



1 + g 2 /ir 



which, again, turns out to be ambiguous. Like we did with the same ambiguities before, let us evaluate the summation 
over n in ( p8[ ) (which is reminiscent of imaginary-time integrations) before the spatial summation and let L — > oo. 
Then one gets 



l+g 2 /ir 



{x x - yi)7 i 



(89) 



VI. CONCLUSIONS AND DISCUSSION 



The main new results obtained in this work are the following: 

1) In the imaginary time formalism, the fcrmionic generating functional Zf with an external electromagnetic field 
and the full generating functional Z for the Schwinger model have been explicitly obtained for any spatial length L 
in the trivial sector (in which the Dirac operator has no zero modes). 

2) The work previously done in at finite T but fi = 0, in which the model was formulated in a two-dimensional 
torus for an arbitrary number of zero modes, can be extended when both T and fi are nonzero. Such an extension has 
to be worked out carefully due to some non trivial peculiarities of the fi ^ case. Technically, the main distinctive 
feature is the lack of hermiticity of the Dirac operator. This implies a non- vanishing phase factor for the fermion 
determinant in the sector with k zero modes. Using functional methods we have evaluated this term for k — (the 
trivial sector) , which plays a crucial role in the solutions for the Schwinger and Thirring models presented here and in 
the physical features thereof. That phase depends on T and fi, is linear in the zeroth component of the electromagnetic 
potential Aq (in agreement with charge conjugation symmetry arguments) and vanishes if fi = for any T, Aq. 
Furthermore, this term is topological, in the sense that it changes only under nontrivial gauge transformations, with 
nonzero winding number around the S 1 parametrising the Euclidean time. In terms of the Hodge decomposition of 
the gauge field in the torus, it only depends on the harmonic part. The existence of topological /i-induced effective 
actions seems to be a common feature of different models [ p3|p4| ] . 

3) For the Schwinger model we have calculated the thermodynamical partition function. The topological phase 
factor in the effective action gives rise to a nonperturbative contribution that in the L — ► oo limit exactly cancels the fi 
dependence contained in the free fermionic partition function. Then, at L — ► oo, the partition function is independent 
on fi and hence the total charge density of the system is zero. In other words, the system bosonizes even though it 
could have a net fermionic charge density at nonzero fi. The partition function factorises, as in || (into that of free 
fermions at fi — times a factor, which is the ratio of that of massive bosons, divided by that of free massless bosons), 
the mass of the boson being m = e/^/n, independent of T and fi. The exact boson propagator has an additional 
/i-dependent piece for any L. At L — > oo this new piece gives rise to a vanishing inverse correlation length squared 
M 2 , which is interpreted through the relationship between M 2 and the first derivative of the charge density given 
in [jl6[ . However, by calculating for /i ^ the thermal average of the Polyakov loop (which is /i-dependent) and its 
correlator (/^-independent), we have shown that the Z symmetry is broken for any T and fi (deconfinement) and that 
the screening mass between two opposite charges is the mass m of the boson. A study of what happens regarding the 
above mentioned cancellation in the thermodynamical partition function, when L is kept finite, lies outside the scope 
of this work. Our computations of the thermodynamical partition function through two different methods yielding 
the same result establish the consistency of our approach. 

4) Several important features of the solution of the Schwinger model for fi ^ in the sectors with zero modes 
are summarised, as they are closely related to the analysis in the trivial sector. Namely, we have given the general 
structure of the fermion determinant, solving the spectrum of the Dirac operator for an instanton configuration when 
fi 7^ 0. The chemical potential breaks the chiral degeneracy of the spectrum. We remark that the correlation functions 
for fi ^ have been analysed in || , although in that paper a different approach based on bosonization is used and the 
harmonic part of the gauge field (and hence the contribution of the phase factors) is not considered. The analysis of 
the phase factors and the fermionic two point function when there are zero modes lies beyond the scope of this work. 

5) In the imaginary time formalism as well, the generating functional for the massless Thirring model at finite T, 
fi is constructed in terms of the fermionic generating functional Zf for the Schwinger model, previously found in this 
work. We have justified that it is enough to restrict ourselves to the trivial sector for Zf. The thermodynamical 
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partition function, the total fermion number density and the fermion correlation function have been computed for 
non-vanishing fi and T. A distinctive feature is that all of them depend nontrivially on /x, as a consequence of the 
non-trivial phase of the Schwinger model. Our result for the pressure agrees with (m) , which shows that our 
different approach is consistent. Our total fermion density differs from jl3|j , where it was obtained, using real time 
formalism, that the model is equivalent to that of free fermions. The origin of that discrepancy is that in [ fi"3| the 
harmonic pieces of the vector field are not considered. 
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APPENDIX A: THE FERMION FUNCTIONAL DETERMINANT FOR $ / 



We shall outline here the calculation of det' H(A; fi) in section II C , following steps similar to those in |Q, with 
suitable generalisations for our present case. Firstly, we shall relate det' H(A; fi) with det' H (A; fj,). For that purpose, 
we define p a replacing e — > ea in (|l7|). Then, the corresponding H a = p a fl a interpolates between H(A) and H(A) 
when a varies from to 1, and so on for H a = TftJP^. The operator H a can be cast as in (|l^), but now with 
cte^ u d v (j) and the operator H a is obtained from H a by changing /i — > —fj,. By using £- regularization [p7| , 



A^ — A^ 
we have: 



log det' H Q = - — ( H {s;a) 
as 



s=0 



Cff(s;a) 



q=k+l 



(a) 



(Al) 



H q (a) denoting, generically, the non- vanishing eigenvalues of H a . As in section II A, we choose the eigenstates of 
H a as ip q a ^ = (l/) a 4> q a ^)l \J 'fi q (a), where 4> q a ^ are the eigenstates of H a for fj, q (a) ^ 0. Now we use the Feynman- 



Hcllmann formula u a (a) = (4> q a \ H a (j) q a ')^ where the dot indicates derivation with respect to a, and the Seele y-de 



do.) 



Witt expansion (27) for H a . Then, following similar steps as in we can write the derivative of logdet'7f Q in (Al) 
with respect to a, in terms of E a = E + aAcjj, 4>{x) and the zero modes (f>^ and tp^ of H a and H a . The latter are 
related to the zero modes of H and H simply by multiplying by exp(— eajs). Then, the integral in a can be done 
and we obtain 

f27t]. p7| for through: 



det'H(A; //) = det'if (A; fx) det <j 7V (1) 
2e 2 



x exp ■ 



d xcj){x) 



E- 



N (o) 
1 



det <^ M (1) 



A(j){x) 



(A2) 



where the elements of the matrices and are N^) = J T cPx/fi^^ 4>^ and M^) = J T d 2 



(")t ; » 



Secondly, we are going to calculate the spectrum of H(A; /x) in (|28[), with the boundary conditions discussed in the 
text. Like in the /i = case 0, we shall try eigenfunctions bearing the form: 



X(fp (p. 



(A3) 



By plugging (A3) into the eigenvalue equation, we arrive at an harmonic oscillator eigenvalue problem, which can 
be solved in the standard fashion. However, the functions (jr 11 in (A3) do not satisfy the right boundary conditions. 
It is not difficult to see that: 
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(A4) 



with k' = n + — n_ = $/27r, are the co rrec t eigenfunctions, which do satisfy the right boundary conditions in (0). We 
have used that the functions in (A3) depend on x 1 and n only through the combination y = x 1 + 2nL(n + 1/2 — 

efto)/*. 

The usual harmonic oscillator quantisation condition for the & m states reads in this case: 



X ± L/3 



± sgn$ = 2m + 1 



(A5) 



m being an integer, with to > 0. Then, the eigenvalues for $ 7^ are independent on /i, and they are: A = with 
degeneracy k, and A m = 2to|$|/L/3 with degeneracy 2k. As we had anticipated, the zero modes appear with only one 
chirality, equal to the sign of <E>. The eigenfunctions are those in ( A4) and (A3) with 



Hn 







[ 1*1 2I 


V l p 


exp 


2L(3 V 



(A6) 



where H m are the Hermite polynomials. Once we know the eigenvalues, we can calculate the determinant using again 
(■-function regularization and we get: 



det' H(A; n) = exp 



as 



«=o. 



= (* L P\ 
\ 1*1 ) 



(A7) 



independent on /i. If we concentrate only in the zero modes, that is, the (n, to) eigenstates in (A4) with to = 0, it 
turns out that they are already orthogonal and that their norms are independent on /z: 



ll^,o|| 2 = 



itL(3 3 



(A8) 



As it was commented in the text, once we know the spectrum of the Dirac operator for $ 7^ we could calculate the 
chiral condensates, which do not depend on the Green function G(x, y, eA; fj,), up to the phase factor. In particular, 
we have 



{il>(x)P±1>(x)) f =exp [iJ w (A ;t i)\ ^det'H(A ;f i)<p\(x)P ± Mx) 
where the superscript 1 indicates the sector with only one zero mode. 
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